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Abstract 



The generalized relativistic effective core potential (GRECP) approach is 
employed in the framework of multireference single- and double-excitation 
configuration interaction (MRD-CI) method to calculate the spin-orbit (SO) 
splitting in the ^P° ground state of the Tl atom and spectroscopic constants 
for the 0"*" ground state of TIH. The 21-electron GRECP for Tl is used and 
the outer core 5s and 5p pseudospinors are frozen with the help of the level 
shift technique. The spin-orbit selection scheme with respect to relativistic 
multireference states and the corresponding code are developed and applied in 
the calculations. In this procedure both correlation and spin-orbit interactions 
are taken into account. A [4,4,4,3,2] basis set is optimized for the Tl atom 
and employed in the TIH calculations. Very good agreement is found for the 
equilibrium distance, vibrational frequency, and dissociation energy of the TIH 
ground state {Re = 1.870 A, lo^ = 1420 cm^^, = 2.049 eV) as compared 
with the experimental data {Re = 1.868 A, Ug = 1391 cm^^, De = 2.06 eV). 
SHORT NAME: GRECP/MRD-CI calculations on Tl and TIH 
KEYWORDS FOR INDEXING: Relativistic Effective Core Potential 
Configuration Interaction, Molecule with heavy atoms, Electronic structure 
calculation. 

31.15.-hq, 31.20.Di, 71A0.+X 
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I. INTRODUCTION 



During the last few years a large number of publications have dealt with calculations of 
the ^-P°/2 ^-^3/2 splitting in the ground state of the Tl atom and spectroscopic constants for 
the 0+ ground state of TIH. Such interest to these systems arises because of their relatively 
simple electronic structure in the valence region. This makes them very convenient objects 
for testing methods for the description of relativistic and correlation effects. We can mention 
some recent papers in which the electronic structure of thallium was studied and papers 
UsHTlj in which the calculation of spectroscopic constants for TIH was carried out. With the 
exception of the atomic RCC calculation by Eliav et al. and the atomic CI/MBPT2 
calculation by Dzuba et al. , the published results cannot be considered to be very accurate 
and reliable, however, primarily because of the rather small basis sets and the small numbers 
of correlated electrons. 

In calculations of Tl and TIH with the use of the relativistic effective core potential 
(RECP) approximation [|1^], in which only 13 thallium electrons are treated explicitly (13e- 
RECPs), one more problem appears. The correlation of the outer core (OC) and valence 
(V) electrons, occupying the bd and ns,np,nd (n = 6, 7, . . .) orbitals, respectively, cannot 
be satisfactorily described, mainly because the smoothed V-pseudoorbitals (pseudospinors) 
have the wrong behaviour in the OC region. One-electron functions (p'^^'k' (r) , being some 
linear combinations of virtual orbitals, correlate to occupied orbitals (p"!^'^ (where x = c,v 
stands for the OC and V orbital indices) and are usually localized in the same space region 
as (j)"'^'^. Therefore, the original "direct" Coulomb two-electron integrals describing the OC-V 
correlation of (j)"'"^ and (j)"'"^ can be well reproduced by those with the pseudoorbitals, despite 
their localization in different space regions. However, a two-electron integral describing the 
"exchange" part of the OC-V correlation. 



» / rfr'Cr^(rOCV)^, (1) 
Jr' r — r 



cannot be well reproduced because the V-pseudoorbitals are smoothed in the OC region 
where the OC-pseudoorbitals are localized (for more theoretical details, see Ref. ||T3|). 

The first RECPs for Tl with the 5s, 5p shells treated explicitly (21e-RECPs) for which this 
disadvantage of the earlier "semicore" RECPs was overcome were generated and tested in 



single-configurational calculations by Mosyagin et al. [0,0 • Some other inherent problems 
of the "nodeless" RECPs were also solved with the 21-electron Generalized RECP (21e- 
GRECP) version presented in Ref. [0,0. In Ref. [|T§, for the case of the 21e-GRECP it 



was also shown that the 5s, 5p pseudospinors could be frozen while still providing significantly 
higher accuracy than 13e-RECPs because the valence and virtual ns and np {n = 6,7, . . .) 
pseudoorbitals in the former case already have the proper nodal structure in the OC region. 



II. THE GRECP OPERATOR IN THE SPIN-ORBIT REPRESENTATION 

In most existing quantum-chemical codes for molecular calculations with RECPs (as well 
as in the MRD-CI code used in the present work) spin-orbit basis sets are used. In these 
versions the number of the two-electron integrals is substantially smaller than in the case 
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of spinor basis sets providing the same level of correlation treatment. Spin-orbit basis sets 
are preferable in the calculations in which correlation effects give a higher contribution to 
the properties of interest than those of a relativistic nature. This is usually the case for 
valence and outermost core electrons, which mainly determine chemical and spectroscopic 
properties of molecules. 

Together with the spin-orbit basis set, the GRECP for Tl should also be employed in 
the spin-orbit representation. Following Ref. [ITBIO], the components of the spin-averaged 
part of the GRECP operator called the averaged relativistic effective potentials (AREP) are 



written in the form 
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where Uni±{r) are the potentials generated for the (pni±ij') pseudospinors by means of the 
Goddard scheme [T^; Ux is the principal quantum number of an outercore (n^), valence (n^,) 
or virtual (n^) pseudospinor; / and j are angular and total electron momenta; ± stands for 
j = / ± 1/2; P„j±(r) is the radial projector on the OC pseudospinors: 



(5) 



Clearly, the AREP component of the GRECP may be used in calculations with nonrela- 
tivistic quantum-chemical codes in order to take account of spin-independent relativistic 
effects. 

The operator of the effective spin-orbit interaction can be derived following the expression 



for the spin-angular projector P;-|- from Ref. [|17 
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Its components, called the effective spin-orbit potentials (ESOP), can be written as |1T3tp^ 

AUUr) = U^j+{r)-Unj-ir), (7) 
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I 

Pi= \lm){lm\ , (10) 



-I 



where \lm){lm\ is the projector on the spherical function Yim- 

Neglecting the difference between U^Jf'^ and t/n„Lj for virtual pseudospinors with I > L 
(for theoretical details see Ref. |13]), one can write the GRECP operator U as 

1=0 

+ EEuir^WP. + E[u„T^ + Euf.f • (11) 

rtc 1=0 1=1 

Note that the nonlocal terms with the projectors on the most important correlation func- 
tions 'f^j±-nki^±{r) (P (where x = c,v) localized mainly in the OC and V regions and with 
the corresponding potentials U^°Ji'^j_{r) can be taken into account in the considered expres- 
sions for the GRECP operator additionally to those with the OC projectors. Obviously, the 
non-local GRECP terms for the frozen OC pseudospinors can be omitted in the sum over 
(nj) in Eq. ([Tl|). 

We should emphasize that in spite of the rather complicated form of the above GRECP 
operator, the main computational effort in calculating matrix elements with the GRECP is 
caused by the standard radially-local operator, which is also a part of conventional RECP 
operators, and not by the non-local GRECP terms. Thus, the additional complications in 
calculations with GRECPs are negligible in comparison with treatments employing conven- 
tional semi-local RECPs if comparable gaussian expansions are used for the partial poten- 
tials. The more critical point is that the effort in the calculation and transformation of 
two-electron integrals is always substantially higher than that in the computation of RECP 
integrals for all known RECP versions (including GRECPs) when appropriately large basis 
sets are employed in the precise calculations. 



III. FROZEN-CORE APPROXIMATION FOR THE OUTER-CORE SHELLS 

To perform precise calculations of chemical and spectroscopic properties, correlations 
should be taken into account not only within the valence regions of heavy atoms and heavy- 
atom molecules but in the core regions and between the valence and core electrons as well. 
In practice, the goal is to achieve a given level of accuracy by correlating as small a number 
of electrons as possible, thus reducing the computational effort. However, as discussed in 
the Introduction, the accuracy of the RECPs generated for a given number of explicitly 
treated electrons cannot always satisfy the accuracy requirements expected from correlating 
all these electrons in the corresponding all-electron calculation. This is true, in particular, 
for calculations of Tl, having a 5d^^6s'^6p^ leading configuration in the ground state, and its 
compounds. 

To attain an accuracy level of 400 cm"^ for the ^-P{/2 ~ ^^3/2 splitting in the ground 
state and for excitation energies to low-lying states of Tl and to take account of the core 
polarization, one should correlate at least 13 electrons, i.e. include the 5d shell. This 
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is achieved in the present MRD-CI calculations with / and g basis functions describing 
mainly polarization of the 5d shell (for other recent results see, e.g., Some data 

from our 13e-CI calculations of the SO-splitting in the ground state of Tl are collected in 
Table | in comparison with the 3e-CI results, which in our DF/CI (Dirac-Fock calculations 
followed by CI) and the GRECP/CI calculations have errors of about 600 cm^^. 

We also should mention the recent relativistic coupled-cluster (RCC) results of Landau 
et al. 1]^, in which 35 electrons are correlated and a decrease of close to 90 cm~^ in the above 
mentioned SO splitting is due to the Breit interaction. Note that this interaction is not yet 
taken into account in the RECPs considered in the present work. 

Obviously, the 5d shell should also be explicitly treated in calculations of molecules 
containing Tl to take into account core relaxation and polarization effects with satisfactory 
accuracy. For these calculations it would be optimal to use the RECPs with 13 electrons 
of Tl treated explicitly (13e-RECPs) such as the RECP of Ross et al. |19| or our valence 



RECP version [|T^ . None of the known nodeless 13e- RECPs can provide the aforementioned 
accuracy, however. Although single-configurational tests [0,0 give errors of 100 cm~^ or 
somewhat more for excitation energies to low-lying states, they are dramatically increased for 
13e-RECPs if all 13 electrons are correlated. The reasons are discussed in the Introduction 
(one can also see the results of the 13e-RECP/MRD-CI calculations in Ref. and of the 
13e-PP/MRCI calculations in Ref. 0]). 

To overcome this disadvantage, one should use RECPs with at least 21 electrons, e.g. 



21e-GRECP ||1^JT5[| and 21e-PP |^ for Tl. The 5s and 5p pseudospinors can be treated 
as frozen, however, while still providing the aforementioned accuracy. The bp orbitals have 
energies about four times higher and their average radii are 1.4 times shorter than those 
for the bd orbitals. Moreover, their angular correlation is supressed as compared with the 
bd shell because the most important polarization functions {bd for the 5p orbitals and 5p 
for the 5s orbitals) are completely occupied in the lowest-lying states. Therefore, the 5s, 5p 
orbitals are substantially less active in chemical processes. 

In order to freeze the 5s and bp pseudospinors, one can apply the energy level shift 
technique 0]. Following Huzinaga et al. [^|, one should add the matrix elements of the 



SCF field operators (the Coulomb and spin-dependent exchange terms) over these OC pseu- 
dospinors to the one-electron part of the Hamiltonian together with the level shift terms 

Y: i?n.,«± Pn.,;±(r) , (12) 

where Bn^^i± is at least of order |2£„^^;±| and en^^i± is the orbital energy of the OC pseu- 
dospinor 0nc^.i±(^) to be frozen. Such nonlocal terms are needed in order prevent collapse 
of the molecular orbitals to the frozen states (the 5si/2, 5pi/2,3/2 pseudospinors for Tl). All 
terms with the frozen core pseudospinors described here (the Coulomb and exchange inter- 
actions, and the level shift operator) can easily be presented in spin-orbit form with the help 
of eq. (H), as was done above for the GRECP operator. 

More importantly, these OC pseudosj»mors can be frozen in calculations with spin-orbit 
basis sets and they can already be frozen at the stage of calculation of the one-electron 



matrix elements of the Hamiltonian, as implemented in the MOLGEP code [^. Thus, 
any integrals with indices of the frozen spinors are completely excluded after the integral 
calculation step. 
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In single-configurational calculations with the numerical HFJ code |T5| we have seen 
that the SO splitting of the 5p shell increases the resulting SO splitting of the ground 
state by about 400 cm~^, whereas the SO splitting of the 5d shell decreases the final SO 
splitting by almost the same value. Therefore, it is important to freeze the 5pi/2 and 5ps/2 
(pseudo)spinors and not some averaged 5p (pseudo)orbitals if the SO interaction is to be 
taken into account in the 5d and valence shells. 

In Ref. 0], the 2 le- "energy-adjusted" Pseudopotential (PP) having the features which 
have been emphasized ||13|- |T5tp2| as inherent for GRECPs (different potentials for the 5p 
and 6p pseudospinors in the case of Tl) is generated and applied to the calculation of the SO 
splitting in Tl, with the core correlations described by the core polarization potential (GPP). 
Some average OC pseudoorbitals are frozen and the SO splitting of 7810 cm~^ obtained in 
their 21e-PP/MRCI calculation is quite different than our result. 

After applying the projection operator of eq. @ to the level shift ([T^) , Coulomb and ex- 
change terms with the frozen core pseudospinors, the AREP and ESOP parts of the GRECP 
operator are to be modified to include these new contributions. This technique was success- 
fully employed in our earlier calculations of the spin-rotational Hamiltonian parameters in 



the BaF and YbF molecules 231. 



The freezing technique discussed above can be efficiently applied to those OG shells for 
which the spin-orbit interaction is clearly more important than the correlation and relaxation 
effects. If the latter effects are neglected entirely or taken into account within "correlated" 
GREGP versions [0], the corresponding OG pseudospinors can be frozen and the spin-orbit 
basis sets can be successfully used for other explicitly treated shells. This is true for the 
5pi/2,3/2 subshells in Tl, contrary to the case of the 5^3/2,5/2 subshells. Freezing the OG 
pseudospinors allows one to optimize an atomic basis set only for the orbitals which are 
varied or correlated in subsequent calculations, thus avoiding the basis set optimization 
for the frozen states and reducing the number of the calculated and stored two-electron 
integrals. Otherwise, if the 5p shell should be correlated explicitly, a spinor basis set can be 
more appropriate than the spin-orbit one. 



IV. THE MRD-CI METHOD 



In the multireference single- and double-excitation GI approach |^ , the AS-basis sets of 
many-electron spin-adapted (and space symmetry-adapted) functions (SAFs) are employed. 
This method makes use of configuration selection and perturbative energy extrapolation 
techniques and employs the Table GI algorithm ||25| for efficient handling of the various 
open-shell cases which arise in the Hamiltonian matrix elements. Some new features of the 
selection scheme used in this work are considered below. The higher excitations in the GI 



treatment has been assessed by applying the generalized multireference analogue ||2^ of the 
Davidson correction |]26| to the extrapolated T=0 energies of each root. 

After selecting the AS'-sets of SAFs for a chosen threshold Tj {i = 1, 2), they are collected 
together in accord with the relativistic double-group symmetry requirements and a spin-orbit 
GI (SO-GI) calculation is performed with these SAFs to obtain some SO-roots (^J^'^') 
and their energies [Sj''') which are of interest in a considered double group irreducible 
representation (irrep). Then the linear T=0 correction is evaluated in the basis of the 
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calculations with the Ti and T2 thresholds. Finally, the generalized Davidson (or full CI) 
correction is applied to each root of interest. 

The stage of the molecular spectroscopic constants calculation begins with the fitting 
of the relativistic CI potential curves to polynomials which are employed to construct ap- 
propriate Born-Oppenheimer nuclear motion Schrodinger equations solved by the Dunham 
method with the help of the DUNHAM-SPECTR code of Mitin \M. 



A. Features of the spin-orbit selection procedure 

Let us define a Hamiltonian H for a molecule as 

H = H(°) + V™"' + H^^, (13) 

where H'-'^-' is an unperturbed spin-independent Hamiltonian, V^°''^ is a two-electron oper- 
ator describing correlations, and H'^*^ is a one-electron spin-orbit operator (ESOP in our 
case). Let us choose an orthonormal basis set of SAFs in the AS'-coupling scheme 

(or "spin-orbit" basis set). In particular, these SAFs can be solutions of Hartree-Fock equa- 
tions with a spin-averaged RECP for the molecule considered. The H*^'') Hamiltonian is 
constructed to be diagonal in the given many-electron basis set: 

H(0)$(")^^ = ^(n)A5^(n)A5 ^ ^^^^ 

where n = 0, 1, . . . (see below the description of the indices in more detail). Additionally 
define H*^°^ so that 

< $(")^^|H(°)|$i")^^ > = < <|.(")^^|H|<|.i'^)^^ > (15) 

in order to exclude the first-order PT contributions to total energies of molecular states (this 
corresponds to the Epstein-Nesbet PT form). 

We will ignore the two-electron spin-dependent (Breit) interactions which ordinarily can 
be neglected when studying chemical and spectroscopic properties. Breit and other quantum 
electro dynamic (QED) effects are relatively large for lanthanides and actinides, but for the 
V and OC shells they can be efficiently represented by the one-electron j-dependent RECP 
terms. 

Let us distinguish the following types of many-electron functions which are considered 
in a double-group symmetry: 

• ^^^^^^^ ^ Ej^^^^ }f^Q^^ are reference SAFs ("Mains") and their energies 

at n = for those AS'-irreps which are of interest for the final spin-orbit CI (SO-CI) 
calculation; 

• ^-q/^p^^ ^gj^^^^yY^Q^'^ are some of the CI solutions ( "AS-roots" ) and their energies 

g{0)^S ^(0)A5|jj(0) ^ YCorr|^(0)A5 ^ ^^^^ 

in the AS'-irrep which diagonalize the (H^^-' + V™^'~) in the subspace of Mains only; 
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• {'^f^^'^ , £^jo)SO|^(o)so gome of the SO-CI solutions ("SO-roots" which are of inter- 
est) and their energies 

which diagonahze the complete H Hamiltonian in the subspace of all Mains collected 
from all the AS'-irreps considered; 

• {^f^^^,Ej^^^^}f^Q^^ are the singly-excited SAFs (SE-SAFs) and their energies (|TB|) 
at n = 1, i.e. 

where P^"^ = [AS >< A^l is a projector on the subspace of the AS'-states, (a^) 
are the creation (annihilation) operators of one-electron states (spin-orbitals) 0p(0q). 
The SE-SAFs can be automatically selected because of their relatively small number; 

• {$(2)A^?^^(2)A5|^j2)As doubly-excited SAFs (DE-SAFs) 

$r ' e {P^%+a>. a. \ U {$rn) V (p, g, r, s; J, K, L) (20) 

and their energies (|T^ at n = 2; a SAF $} should be selected in accordance with 
some selection criteria to be used in the final SO-CI calculation. In principle, triple 
and higher excited sets of SAFs can be similarly defined. 

The correlation operator, V™^'', has the symmetry of the molecule and, therefore, can 
be rewritten as 

-^corr — "y ^ pAS-^corrpAS (21) 
AS 

It normally gives the most important contribution through the second-order Brillouin- 
Wigner PT energy correction in the basis set of (after appropriate redefinition of 

H^°) in the subspace of Mains, see Ref. [p^ ,^): 



I ^ fT,('")A5'|ArcorrhT,(0)A5 ,2 

E>p I ^ I ^ 1^0 ^1 

^ cAS rp{n)^S 

n=l,2 J Cn — i^T 



(22) 



for the non-degenerate ground state \E'q'^ with the exact energy in the AS'-irrep (ob- 
viously, terms with n > 3 are automatically equal to zero because V™*"^ is a two-electron 
operator). A similar expression with the replacements \E'q°'''^'^ — "^f^^^ and S^^ — > S^^ can 
be applied for excited states (some precautions should be taken concerning the degen- 
erate states and the orthogonality constraints with respect to the lower-lying states with 
J < I). As a result, the first rows, columns and energies on the diagonal of the Hamiltonian 
matrix 

< $W^^|V--|#^^ > , < $J^)^^|H|<|.!r^^^ > (23) 
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for n = 1,2 are usually employed in the selection procedures for SAFs based on 

the nonrelativistic and approximations (when H"^*^ is not taken into account) P5|J5D | 



or on the multi-diagonalization scheme for subsequent calculations of "^f^ . In spite of 
some differences between these selection schemes, they are not very essential for the final 
CI results if a high quality reference set (set of Mains) and a suitably small threshold are 
chosen. 

For molecules with heavy and very heavy atoms, the H'^*^ operator can give large contri- 
butions to the energy both in second and in higher PT orders if a non-optimal set of Mains, 
{<l>j°''^'^}, is chosen after an SCF calculation with the SO-averaged potentials (AREPs). The 
latter is the usual practice and the set of Mains generated in such a manner can be smaller 
than optimal for the case of large SO interaction. Therefore, not only second but third and 
maybe even higher PT order(s) can be important in the selection procedure for a "bad" 
set of the starting roots "^f^^^ . This means that the off-diagonal matrix elements of H be- 
tween secondary many-electron basis functions (SE-, DE-SAFs) may be introduced into the 
selection procedure because H'^'^ is a substantially off-diagonal operator contrary to V™^^: 

< ^("^^^IH^^I^J,"'^^'^' > (A^) and (A'^') can be different, n' e {n, \n ± 1\} . (24) 

In particular, H'^'^ gives zero matrix elements between SAFs belonging to the same AS'-irrep 
in the or 6*2^ symmetry groups. 

For simplicity, let us consider the selection scheme based on the approximation (p2|). 

In the nonrelativistic-type selection scheme, a SAF $ j is selected in a AS'-irrep if 



< $(l,2)AS|^eor-r|^(0)A5 ^ 



> 5E^^ , (25) 



where / < A^'^^ and 5E^^ is a threshold criterion for the energy selection scheme in the AS'- 
irrep. In (p5|) we have replaced the exact energies by the approximate g'f'^^^ values (p!7|) 



that corresponds to the Rayleigh-Schrodinger PT case. Such a simplification is justified for 
small 5E^^ and good reference states. 

In a SO-CI calculation within some relativistic double-group irrep, substitutions for the 
reference state {mf^^ ^J°)^^) and the perturbation (V™^^ ^ V^"'"'^ + H^^) should be 
used in the previous expression, so that 



,(l,2)A5|^^corr _L w5OhT/(0)'S'O ^ |2 

> 5E^^ , (26) 



t;i(1,2)A5 c{0)SO 

where SEj^'^ is a selection threshold for ^^j'"^^^^ to be used in the subsequent SO-CI calcu- 
lation. 

In more detail, the matrix element in the PT numerator of the above formula can be 
rewritten as 

I < $a,2)A5|^eor-r|^(0)SO ^ p (27) 

+ 2K(< vl/f ^«|V— X $«^^|H^O|vl/f ^« >) (29) 
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by taking into account eq. (pTD in the calculation of the matrix elements for V'^°''^, contrary 
to those for H"^*^. In spite of mixing different AS'-states due to H*^*^, the number of non-zero 
matrix elements with H"^*^ in eq. (^) is usually relatively small because the SO interaction 
is a one-electron operator (see eq. (pi]) ) which is very localized compared with the long-range 
Coulomb interaction. Thus, one can see that the nonrelativistic-type selection due to V'^"'''' 
with respect to {P^'^vE'^'''*'^*^} in each considered AS'-irrep and automatic selection of all SE- 
SAFs (H) can be efficiently applied instead of eq. (|26|) . It must be emphasized that 

contrary to the selection schemes in the nonrelativistic case, SE-SAFs should be generated 
with respect to the Mains from all the used A'S"-irreps. In a more simplified treatment, the 
automatic selection of SE-SAFs can be done with respect to a subset of the most important 
Mains, e.g. having largest Cl-coefficients in the '^f^^'^ roots. 

Next let us consider the terms from the third-order PT energy (PT-3) for SAFs '^^^^} 
which can be essential for the SO selection procedure. Below we shall discuss only matrix 
elements in the PT numerators of the corresponding PT-3 terms because specific expressions 
for the energy denominators are not essential for our analysis and conclusions. For simplicity, 
we shall omit the terms conjugate to those considered. 

The first two types of the PT-3 matrix elements are: 

< xi.W^O|H^«|$?)^^ X $«^^|H^O|<|.«^'^' X <|.«^'^'|H^O|^f)^« > , (30) 

The first intermediate state, ^^p^^ , is a test SE-SAF and the indices for other intermediate 
SAFs run over all the allowed ones. The PT-3 terms summed over the indices of the second 
intermediate state give contributions (together with the conjugate terms) for the selection 
of the test SE-SAF. However, the SE-SAFs can be selected automatically and these terms 
are out of our particular interest. 
The following matrix element type 

< ^0|^eorr|^(l,2)A5 $a,2)A5 1 jj50 1 ^(l)A'S' ^^A'S' | jjSO|^(0)SO ^ _ ^32) 

can be used for the selection of $j '^^^'^ and ^^^^ ^ when summing over another set of 
intermediate states in the PT-3 expression. As one can see, this term can be used for the 
selection of both SE-SAFs and DE-SAFs. The above expression is quadratic in the (large) 
H'^*^ interaction contrary to the remaining terms considered below. The contribution of the 
terms with matrix elements (|32D can be essential and their use for the selection of DE-SAFs 



^(^2)A5 important for a subsequent SO-CI calculation. 

The following matrix element types contain a second order perturbation in V™'^'' and, 
therefore, we can suggest that in general they are less important for our consideration than 
the above terms: 

< ^(0)^^O|^,o,,|^{l,2)AS ^{l,2)AS|^eorr|$(l)AS ${1)A5 | jjSO| ^{0)SO ^ ^ ^33^ 

< ^(0)SO|V™-|$(/'2)A5 ^(l,2)A5|jj50|${l,2)A'5' ^{1,2)A'S' |^,or-r | ^(0)50 ^ _ ^3^^ 

These terms, together with the conjugate ones, can be used for the selection of $y and 
^{i)A5^ The term 
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can be analyzed separately because it contains both the intermediate states as DE-SAFs. 
In general, it is more difficult to take such terms into account in the selection procedure, 
because of the large number of tested DE-SAFs. We should note, however, that when a 
tested $j DE-SAF is fixed, the other intermediate states, }' those DE-SAFs 

which are only singly excited with respect to the tested one. Therefore, the number of them 
will not be very high. 

For completeness, the matrix element type which is cubic in the V^°''^ perturbation 
should be listed: 

This term is of nonrelativistic type and it is out of our particular interest because it does 
not contain the H'^'^ perturbation. Again, we can separate the term 

from the previous one only because the latter contains both the DE-SAF intermediate states. 

We should emphasize that the terms containing SE- or DE-SAFs in the intermediate 
states of the PT-3 expressions are not taken into account in the and multi-diagonalization 
selection procedures, although these schemes include, in fact, contributions of higher than 
the second-order PT terms. 

When analyzing the above PT-3 terms, it can be concluded that if one replaces the 
reference SO roots, '^^^^'^ ^ by new reference states, ^ which diagonalize the complete 

Hamiltonian H for the sets of both Mains and SE-SAFs taken together, and applies the 
selection criterion based on the second-order PT (p6D, then the main part of the above PT-3 
terms will be taken into account in such a selection. An exception occurs for terms (^) and 
(^), but in general they are thought to be less important than the other third-order PT 
terms. 

In a more sophisticated treatment, the reference '^f^^ ^^'^ SO states can be generated 
when diagonalizing H for the sets of Mains and those SE-SAFs ({$j }), which are au- 
tomatically generated with respect to the most important subset of Mains ({$} }). The 
latter subset can be selected from a preliminary CI calculation for the set of Mains, e.g. 
in a basis of configurations with the highest CI coefficients in '^f^^^ ^ and so on. This is 
worthwhile in order to reduce the number of SAFs in the resulting reference states ^if^^ ^^'^ 
rather than in \^^^^^^^'^ ^ thus reducing the selection time which can otherwise be very large. 

We should also note that the trial SE- and DE-SAFs, which are tested in the above 
selection procedure, are generated only for the set of Mains and not for the {$/°^^ 
Therefore, the number of the tested configurations and the selection time are reasonably 
limited. If the number of configurations used in {$/°^^ '''^'^} is not high, one can extend the 
set of Mains by including the above subset of SE-SAFs, thus obviously enlarging the set of 
the consequently generated and tested SE- and DE-SAFs. 

Again we should emphasize that it is not necessary to use the third-order PT or the 
suggested automatic selection of SE-SAFs in a selection procedure if a fairly good set of 
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the reference roots "^f^^^ is used, i.e., if they provide good approximations to the required 
solutions ^ff^. In particular, if the {'^f^^'^} set is obtained from a preliminary series of 
SO-CI calculations of the studied states, this can be superfluous. 

As an alternative to the above selection schemes with respect to the PT energy, the PT 
expressions for the CI coefficient of a trial SE- or DE-SAF can be also explored. Applying 
the above PT analysis to the case of the '^f'^^ > reference state, a $j ''^^^^ SAF is 
selected if its CI coefficient Cf ''^^^^ satisfies the inequality 

|^(l",2)A5| ^ 

\^ J \ — ^min 1 



where Cmin is the selection threshold for the CI coefficients and 

,(l",2)A5|jj|^(0+l')5O 
^jO+l')SO _ 



^(1",2)AS _ < "^J > foo\ 

-(0+l')5O j-i(l",2)AS ' V-JOJ 



is the first-order PT value for the CI coefficient of a tested SAF which is not included in the 
subset of the $j -'^^ reference SE-SAFs. 

Such a means of selection can be preferable if those properties of primary interest cannot 
be calculated from potential energy curves or surfaces. Moreover, the PT selection with 
respect to both the energy and the CI coefficients can be applied simultaneously if the 
properties are of different nature. 



V. CALCULATIONS 



In the CI calculations of Tl and TIH we used the MRD-CI package combined with 
the SO selection codes based on the scheme described above. Our test calculations have 
shown that spin-orbit selection is very helpful for preparation of appropriate sets of Mains 
and for reducing effort in the final CI calculations with an optimal set of selected SAFs. 



A. Spin-orbit splitting in the ground state of Tl 

Calculations for the Tl atom were performed to optimize the basis set and the level shift 
GRECP parameters for the 21e/8fs-GRECP, i.e. the 21 electron GRECP with 8 electrons 
occupying the frozen OC pseudospinors, 5si/2 and 5pi/2,3/2- The quality of the generated 
basis set is analyzed by calculating the ^-P°/2 ~ ^-^3/2 splitting for the ground state. 

Before discussing the present results, it is worthwhile to to make some brief comments 
concerning numerous values for the T1(^P°) spin-orbit splitting calculated and published in 
the last years. It is well known that this quantity calculated at the one-configuration Dirac- 
Fock level agrees very well, within 100 cm~^, with the experimental value of 7793 cm~^ 
^|. However, such a good agreement results from the fortuitous cancelation of a number 



of large errors caused by the DF approximation. The situation changes dramatically when 
even the three outermost 6s^6p^ electrons are correlated. Some T1(^P°) splitting values 
from our 3e-CI calculations employing different codes and basis sets are given in Table | 
and lie between 7130 and 7210 cm~^ The corresponding 3e-CI results obtained by other 
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groups after 1996 range from 6800 to 7800 cm~^ and such a large divergency can not be 
considered as satisfactory, because the ground state of Tl has a very simple configuration 
structure as compared to other heavy elements. We consider our calculated values of about 
7200 cm~^ for this splitting as reliable for an approach in which the 5d spinors are frozen 
after the DF calculation of the nonrelativistically averaged 6s^6p^ configuration. Taking 
into account that a contribution of aproximately -100 cm~^ arises from the Breit terms, the 
deviation from the experimental value for the splitting is around 600-700 cm~^ (this size 
of an error can be justified theoretically). As will be shown below, the computated value 
for the T1(^P") splitting can be significantly improved if 5d electrons are explicitly included 
in the calculations and the corresponding basis set contains functions with sufficiently high 
angular momenta. 

The optimal basis set was selected in a series of MRD-CI calculations for Tl (with different 
sets of primitives and numbers of contracted s,p, d, f and g functions) to minimize the sum 
of energies for the ground and ^^3/2 states. In these calculations, the SAFs were 

selected in the ^-B2n, and ^B^u irreps of the group (nonrelativistic-type degenerate 
2p ground states belong to these irreps) because these doublets are strongly mixed by the 
SO interaction, resulting in the splitting of the ground state. We have found that two g 
functions should be added to the basis set, giving a contribution of about 9000 cm~^ to the 
2po g^Q^nd state total energies. The resulting [4,4,4,3,2] basis set and GRECP parameters 
for Tl can be found on \http:/ /www. qchem. pnpi. spb.rt\ . 

For the [4,4,4,3,2] basis set we have also performed MRD-CI calculations including SAFs 
from the "^Au irrep and SAFs with quartet multiplicity (^-Biu, ^i?2u5 ^-83^, and ^A„). In our 
calculations with different basis sets, their contributions have decreased the SO splitting by 
about 170 cm~^ and the total energy by about 2000 cm~^. One can see from Table | that this 
decrease is mainly caused by the p- and (i-components which arise from reexpansion of the 
leading spinor configuration in terms of the spin-orbit configurations. For good accuracy we 
can recommend the inclusion of A' IS* ± l|-irreps for the calculation of states having leading 
configurations in AS'-irreps. 

In Table | some of our final MRD-CI results are collected together with the atomic 
relativistic coupled-cluster (RCC) results ||^ obtained with a very large basis set. In these 
MRD-CI calculations altogether 627 Mains in three basic irreps and about 100 Mains in 
five additional irreps were involved and SE-SAFs were automatically generated for three 
Mains to prepare the reference ^^'^}^=i states. Relatively small thresholds, Ti=0.03 

and T2=0.01 fiEh, are used in the final runs with the [4,4,4,3,2] basis set (for the T=0 
threshold and fuU-CI extrapolations), thus selecting respectively about 190000 and 450000 
SAFs altogether. 

One can see that the best SO splitting calculated in the present work underestimates 
the experimental result by about 400 cm~^ (recall that additionally about 90 cm~^ 
is due to the Breit interaction |0). Analyzing our previous GRECP/RCC calculations of 
Hg [^ it can be concluded that this occurs due to the neglect of the OC-V correlations 
with the OC 5p and 4/ shells, and to a lesser extent with 5s rather than due to the atomic 
basis set incompleteness, the GRECP errors or the restricted CI approximation. The OC- 
V correlation (contribution to the total energy) in Tl and Hg will have the same order of 
magnitude for respective pairs of correlated electrons (spinors). 
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We also studied the reliability of the linear T^O extrapolation procedure currently used 
in the MRD-CI code. In the final results of our MRD-CI calculations the corresponding 
correction gives the highest contribution to the cumulative error. So this is a bottleneck of 
the present Tl and TIH calculations with the large number of Mains. 



B. Spectroscopic constants of the ground state in TIH 

The explicit treatment of 5d electrons in precise TIH (TlX) calculations is necessary not 
only due to the strong correlation between these and the valence electrons of Tl, but also 
because of the substantial infiuence of relaxation-polarization effects in this shell on the bond 
formation. This cannot be very accurately taken into account by employing a polarization 



potential [0,0 in combination with, e.g., 3e-RECPs [|T4|,|l5|. The influence of other atoms 
(X) in a Tlx molecule on the 5p, 5s and 4/ shells of Tl is signiflcantly smaller and can 
be neglected if an accuracy of a few hundreds of wavenumbers for the excitation energies 
of low-lying states is sufficient. We neglected their contributions in calculation of the TIH 
spectroscopic constants. 

In calculating spectroscopic properties for the TIH ground state (Table H) we used the 
contracted [4,4,4,3,2] basis set for thallium discussed above and the [4,3,1] set for hydrogen 
(see \http : / / www . qchem.pnpi. spb.riii ) contracted from the primitive (6,3,1) gaussian basis set 
of Dunning ||3^. The SAFs were selected in the ^Ai, ^Bi, and AS-irreps of the 
group because the triplet states are most strongly admixed by the SO interaction to the 
nonrelativistic ^Ai (or in Coov) ground state producing the relativistic 0+ ground state 
in the double C^^, group. 

We have performed three series of TIH calculations for 16 interatomic distances. In these 
runs, the '^'q^^ ^"^'^ reference SO states are generated with the MRD-CI code by diagonalizing 
H for the set of Mains and the SE-SAFs which are automatically selected with respect to the 
single configuration SCF ground state (calculated with the SO-averaged GRECP), giving a 
contribution of more than 90 % to the final wave function. 

The first run is used for preparing an optimal set of Mains for the second series of SO-CI 
calculations. Only one SCF configuration which has the lowest energy in each AS-irrep is 
included into the subspace of AS* Mains and, consequently, the SO reference state consists of 
these SCF configurations and the automatically selected SE-SAFs with respect to the SCF 
configuration from the ^Ai irrep. 

Those SAFs were selected as Mains for the second run which had the highest CI coeffi- 
cients in the first run. As a result, 37 Mains in all irreps together are employed in the second 
run. Relatively small thresholds, Ti=1.0 and T2=0.1 /i-E/i, are used in the second run (for 
the T— extrapolation [0]), thus causing about 20000 and 85000 SAFs to be selected in 
the AS'-irreps altogether. 

In the most computationally consuming third run (with the set of Mains consisting of the 
SAFs having the largest CI coefficients in the wave function from the second run), about 320 
Mains are used altogether and the thresholds are set at Ti=0.1 and T2=0.05 ^lEh. About 
70000 and 130000 SAFs, respectively, were used in the A5'-irreps altogether in the final 
SO-CI calculations. 
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One can see from Table H that the basis set superposition error (BSSE) (see p9| and 



references) must be taken into account for an accurate computation of spectroscopic con- 
stants. The BSSE was studied in the Tl"*" ion calculations for the same interatomic distances 
as in TIH and estimated also in the Tl~ calculations for three distances, i.e. with the ghost 
H atom. The same molecular basis set as in TIH was used for both the Tl and H atoms. 
The contribution from BSSE to the total energy is decisive for the 5d^^ and shells con- 
sidered in the case of T1+, while its changing due to addition of the 6p electrons (which are 
bonding in TIH) can be considered as relatively small, because the difference in BSSE for 
Tl"*" and Tl~ is not significant in comparison with other errors. In the calculations of the 
spectroscopic properties with the counterpoise corrections (CPC), the calculated TIH points 
on the potential curve were corrected with the calculated BSSE for T1+, i.e. for the 5d,6s 
shells taken into account. 

One can see that after applying the T=0, FCI, and counterpoise corrections, the calcu- 
lated properties are in very good agreement with the experimental data both in the second 
and third runs. The accuracy obtained is notably better than for other existing results for 
TIH (and not only for those presented in Table y). We suggest, however, that the very good 
agreement of the calculated with the experimental value can be fortuitous and the "real" 
(full CI) value can be notably different from the listed one because of the approximations 
made. 



VI. RESUME 

The SO splitting in the ground state of Tl is calculated by the MRD-CI method 
with the 21e-GRECP when 5d^^, 6s^ and Qp^ electrons are correlated and the and 5p^ 
pseudospinors are frozen in the framework of the level shift technique. A [4,4,4,3,2] basis set 
is optimized for Tl and an underestimation of about 400 cm~^ is found for the SO splitting 
as compared with the experimental data. 

Further improvement of the accuracy can be attained when correlations with the outer 
core 4/, 5p and 5s shells of Tl and Breit effects are taken into account. We expect that 
this can be efficiently done in the framework of the "correlated" 21e/8fs-GRECP version 
in which 13 electrons are treated explicitly as in the present calculation. The inclusion of 



h-type functions is also desirable, as has been demonstrated for Hg in Ref. |^ . 

Fourteen electrons are correlated in the calculation of spectroscopic constants for the O"*" 
ground state of TIH and very good agreement with the experimental data is found. 

The developed spin-orbit selection scheme and code are demonstrated to be efficient 
when large sets of basis functions and reference configurations are required in high-precision 
electronic-structure calculations. 
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TABLES 



TABLE L Calculations of the spin-orbit splitting of the ground state in Tl (the [. ..],5d 
spinors are frozen from the SCF calculation of the nonrelativistically averaged [. . .]6s'^6p^ configu- 
ration) . 



Method 




SO splitting in cm ^ 


Spinor basis sets: 


[7,7,5] 


f7 7 5 31 


f7 7 5 3 11 


81e-DF/3e-CI 


7129 


7182 


7206 


21e/18fs-GRECP/3e-CI [Ol 

/ / Lt_jJ 


7133 


7187 


7211 


Spin-orbit basis sets: 


[4,4,4] 


[4,4,4,3] 


[4,4,4,3,2] 


21e/18fs-GRECP/3e-MRD-CI (Full CI) 

(^Biu,^B2u,'^B3u irreps of D2h) 


7305 


7373 


7398 


(. . . + '^Au, ^Biu, ^B2u, ^Bsu, ^Au) 


7133 


7205 


7230 


21e/8fs-GRECP/13e-MRD-CI + T=0 + FCI 
(^Biu,^B2u,'^B3u irreps of D2h) 


7332 


7222 


7517 


(• • • + '^Au, ^Biu, ^i32u, '^-Bsu, '^Au) 


7146 


7044 


7380 


Spinor basis set: 






[35,27,21,15,9,6,4] 


81e-DF/35e-RCC § 






7710 



Experiment |3l| 7793 
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TABLE II. GRECP/MRDCI calculations of the spectroscopic constants for the ground state 
of TIH. 







tie 




n 


IVieTjilOQ 






{cTfi 


{ev } 


CO/^TTPV. '~ri [Q Q Ol 1 XJ f/l Q 1 1 

oUUIiIjA: 11 [o,o,o,zj + rl [4,o,iJ 










(^xiaKowitz cz iviarian, lyy* |ioi|j 




i.oD 


iooD 


z. io 


-1 o T~» 1 — ' / O /"A /^T rill f/i A /ill 1 TTf/frxl 

13e-RECP/S0CI: Tl [4,4,4, IJ + H [4,2J 










(UiLabio & Christiansen, 1998 




1.912 


1341 


1.908 


13e-REP/KRCCSD(T): Tl [4,5,5,1J + H [3,2J 










(Lee et al, 1998 |0|) 




1.910 


1360 


2.02 


21e-REP/KRCCSD(T): Tl [4,5,5,1] + H [3,2] 










(Han et al, 2000 pl|) 




1.877 




2.00 


21e/8fs-GRECP/14e-MRD-CI Tl [4,4,4,3,2] + 


H [4,3,1] 








(Present calculations] 










37 Mains, T=0.1 




1.858 


1481 


2.03 


h UPU 




1.8 ^ 2 


1440 


1 no /I 


+ T=0 + FCI 




1.858 


1453 


2.10 


+ T=0 + FCI + CPC 




1.872 


1410 


2.026 


320 Mains, T=0.05 




1.866 


1408 


2.23 


+ T=0 + FCI 




1.858 


1449 


2.124 


+ T-0 + FCI + CPC 




1.870 


1420 


2.049 


Experiment (Grundstrom & Valberg, 1938 [^ 


1) 


1.866^^ 


1390.7 


2.06 


Experiment (Urban et al, 1989 [||]) 




1.872'' 


1391.3 





^Huber & Herzberg (1979) [^ have published value 1.87 A which can be obtained from the 

rotational constant Be. 

^This value is calculated by us from Be. 
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